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ABSTRACT 

Analyses are given for the dynamic response to 
operational maneuvers for spinning symmetric space- 
craft with radial and axial booms. The research 
was performed as part of the prelaunch dynamic 
analysis of the ISEE-3 spacecraft placed in a halo 
orbit around an Earth-Sun libration point, and 
lateY- renamed ICE when it was directed to fly-by 
comet Giacobini-Zinner. The results presented use 
simple spacecraft models with the aim of developing 
understanding, and they frequently gave predictions 
that were good and easily obtained when the results 
from using a general purpose multi -body dynamics 
program were very time consuming to obtain. The 
operations encountered during the spacecraft history 
that are analyzed here are, deployment of radial 
booms, spin-up after partial deployment, station 
keeping, and trajectory changes. The latter two 
can involve both axial thrusting and pulsed radial 
thrusting once per revolution. 

Keywords: Flexible spacecraft, attitude dynamics, 
deployment, spin-up, vibration modes during thrusting. 


1. INTRODUCTION 

This paper addresses a series of problems associated 
with the dynamic response of flexible spinning space- 
craft during thrust maneuvers. It is a companion 
work to Reference 1, which determined spacecraft 
vibration modes and frequencies, and dynamic response 
during spin-up, of the IMP-0 (Interplanetary Monitor- 
ing Platform) spin stabilized spacecraft with radial 
wire antennas. The work reported here was done as 
part of the pre-launch dynamic analysis of the ISEE 
(International Sun-Earth Explorer) series of space- 
craft (whose name had just been changed from IME at 
the time of the study). This series represented a 
follow-on program to the IMP series, employing very 
similar spacecraft configurations. Particular 
attention is given the ISEE-3 spin stabilized space- 
craft with four long radial wire antennas, and short- 
er spin axis booms. 

The spacecraft was subjected to various thrust 
maneuvers while the radial and axial booms were 
deployed, and this paper aims at developing simple 
models to predict the associated vibrational resp- 
onse of the flexible vehicle. The maneuvers con- 
sidered include the following. After orbit injection 
the radial wire antennas are deployed, and the effect 


of the deployment rate history on spin rate and 
spin rate vibations is studied. In order to main- 
tain at least a minimum spin rate, deployment is 
periodically interupted for a spin-up maneuver. 

After the radial and axial wires and booms are de- 
ployed, orbit maintenance maneuvers are required 
to maintain the satellite in its unstable halo orbit 
about the libration point. Because of reliability 
considerations it is undesirable to retract and 
later extend booms and wires for such maneuvers. 

In addition to the orbit maintenance maneuvers, 
there are orbit change maneuvers to go from the 
libration point orbit to a large Earth orbit, and 
eventually to the comet fly-by after various moon 
fly-bys for gravity assist. The satellite has the 
capability of maneuvering in the plane of the spin 
using pulsed radial thrusting that is turned on for 
a portion of each spin cycle. Maneuvering along 
the spin axis can be accomplished with spin axis 
thrusters, or by performing an attitude change and 
using the radial thruster. 

The aim of the analysis here is to use simple models 
to obtain an understanding of the spacecraft vib- 
rational response to these maneuvers. The pre- 
launch dynamic analysis also included many computer 
simulations using large general purpose multibody 
dynamics programs. Often, such simulations were 
quite time consuming, while the equations derived 
here predicted the results accurately and quickly 
in terms that were easily understood and gave more 
insight. 

2. SPIN-UP AND DEPLOYMENT MANEUVERS 

The spacecraft is modeled as a cylindrical rigid hub 
with radial and axial booms represented as rigid 
spherical pendulums, with a linear spring constant 
restoring torque at their pivot points. The spring 
constant is chosen to match the first bending mode 
frequency of the boom as a cantilever beam. This 
simple model makes many of the problems tractable. 

In one case it was possible to obtain results for 
the first vibrational mode starting with partial 
differential equation models, and this resulted 
in the same ordinary differential equation as the 
above simple model, except that the gain coefficient 
in the driving term was altered from 1.5 to 1.566. 
This small difference gives confidence in the other 
results obtained from this simple model. 

Reference 1 develops the fundamental spacecraft mode 
shapes for the satellite before the spin axis booms 
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are deployed. Both deployment and spin-up directly 
excite only the oc mode shown in Fig. 1, in which 
all radial booms simultaneously lag or lead, while 
the hub rotates in the opposite direction. Only this 
mode will be considered here. 

2.1 Derivation of differential equations 

This section derives differential equations in suf- 
ficient generality that they can subsequently be 
used for both spin-up and for deployment, in the 
former case, there is a torque applied to the space- 
craft hub, and in the latter there is kinetic energy 
associated with the radial motion of the booms. 
Figure 2 gives the nomenclature used, where x,y are 
inertial, 0 is the angle of rotation of the hub 
and q is the radial position of a generic volume 
element on a boom of lenght J(t). The portion 
of each boom remaining in the hub at any time is 
modeled as a point mass at the radius r of the 
hub. The total kinetic energy is the sum of: 

1) The kinetic energy of the hub l/2I H Q.6fand 
that of the coiled booms. 2) The kinetic energy 
of the tip masses of mass iil. 3}The kinetic 
energy of the booms. The potential energy of 
all radial booms together is 1/2 koc 1 . Let/° 
be the density per unit length of the booms 
and let their total length be Jt . Then by 
Lagranges equations one obtains the oc equation 

I = I M + M,r* + Z^rL + M-ji* 

M 2= %pL) 

oH) 

[tAjivn^rcoioc] 1 6 + +• [tyX s(« <*-} e~ 

The © equation has a first integral, Eq. 3, when 
no torque is applied to the hub. C is a constant 
of integration, here evaluated for zero initial con- 
ditions. The Lagrange equation for © is Eq. 4 
below. 

[l-an 3 .ri(l- c « , <.“OJe + [.nji + MjrcoicZjjti ( 3 ) 

+ fA| rX>otnoC = C 

C = 1(0) <J 0 = (I Ho + Hf +• *< r "t r z ) 

[x-itV-JUi-caooOle + 

+ C 2 fAirjlc« «]oc +£x- 2 tyi(i- CMoO-aiylWoxJe 
+ Wot _ _ 

X = z t»v + n 2 i]i 


Here Q is the spin-up torque, and u) 0 the initial 
spin rate. 


2.2 Radial boom dynamics during deployment 

We are interested in isolating a differential equa- 
tion for radial boom deformation oc alone, by elim- 
inating 0 dependence in Eq. 2. We assume the boom 
deformations are small so that linearization inoc 
and °c about zero is appropriate. There is no need 
to pick a nominal 0 history to linearize about. 
Decoupling of the linearized equations is complicated 
by the presence of time dependent coefficients, so 


that methods such as Laplace transforms cannot be 
used. Making use of all three equations, Eqs. 2-4, 
one can obtain 

I tCMji 2 - + [RlVlJl -Ir-l 2s (.3M,c*4n 2 J.) 

+ 2(n 1 r + m a 4)Mlj]i<i. + Ll'A + fVic’; 

- C 2 .IM 1 M 2 rX. 4 - HiM 23 pr r)i & 

-XM,H 23 rjL]0C 

= -2 LrM 2 t- ( ia,i r+t\i)n i2 }ic cs) 


This equation is too complex for analytic solution. 

If one considers a constant deployment rate, then 
the coefficients become polynomials in t, but the 
degree of the polynomials is„high. For example, the 
coefficient of 6 L contains t . 

If the deployment rate is slow one could conceive of 
a quasi -static steady state associated with the par- 
ticular solution. Then the approximate particular 
solution is given by dividing the right hand side 
of Eq. 5 by the coefficient of oc on the left hand 
side. This solution often gives rather good results. 
It can be simplified by making further assumptions 
in order to obtain the offset angle during spin-up 
that was derived by approximate means in Ref. 2: 
oa p = - ZL/ (r« 0 y 

2.3 Dynamic response of radial booms during spin-up 

To obtain a differential equation for the radial 
boom deflection oc during a constant torque Q spin- 
up, linearize Eqs. 2 and 4 for oc about zero, and 
for e about the nominal spin rate that would exist 
if the spacecraft were rigid. Decoupling the 
equations, introducing a new independent variable 
T= I +e-t and denoting differentiation with respect 
to x by a prime results in the following equation 
for oc 

oc" + (A.+ r*-)XI^oc = - K (6) 

t = a/(iu.) 

h - -k/(tt 2 .rX.w,j L ) 

fll - { I A o!-)/t6 (I n 3 X 2 - M 2 )] 

K = QMjj/Lefrtvt'-ni)] 

The quantity 6 is small for the spacecraft con^. 
sidered, since the spin-up torque Q is small. SL^ 
is the radian natural frequency of the oC spacecraft 
mode when spinning at cj„ . 

In the companion paper. Ref. 1, Eq. 6 was solved 
analytically for the case of wire booms with k=0. 

The contribution here is to study the effect of a 
nonzero boom stiffness on the solution to Eq. 6. 

In Ref. 1 a transformation was found to convert the 
equation to a Bessel equation, for which the comp- 
lementary function is in terms of Bessel functions 
of order ±1/4. The particular solution for a con- 
stant torque Q was found to be expressible in terms 
of a Struve function of order -1/4. An alternative 
expression for the particular solution was given 
as an infinite series in Bessel functions. Simpler 
approximate solutions were also obtained, by the 
WKB method (Ref. 3) for the complementary function, 
and in the form of an asymptotic expansion for the 
particular solution. 
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It is found hens that the introduction of the boom 
stiffness term k precludes obtaining analytic sol- 
tions in terms of Bessel and Struve functions. 
However, the WKB method and the asymptotic expan- 
sion method can still be applied, although the 
results are more complicated. 

2.3.1 The MKB approximation to the complementary 
function " "" 


Generating the WKB approximation as in Refs. 1,3, 
one obtains the complementary function, Eq. 7 

M. + (felt 1 / 1 )] 4 & } (7 ^ 

>5- z _ Tf Cs) 

whiGh is a good approximation whenever Eq. 8 is 
satisfied which it is in this case because of the 
€* divisor in JT.^. A and i are arbitrary cons- 
tants. It is of interest to determine whether this 
approximate solution is sufficiently accurate that 
there is no large accumulation of phase error after 
many oscillations. Let C + be the arbitrary const- 
tants associated with the'al ternative exponential 
form of Eq. 7. The error in the approximation can 
be characterized by these constants becoming funct- 
ions of time r. An approximation to the total 
change AC in C going from 7:=1 to oo gives 
(Refs. 1,3') 

|AC*/C t \ 6 7/C8il«) (l) 

The bound could be made much tighter, but the pre- 
sence of € as a divisor in nevertheless makes 
this bound small. We conclude that the approximate 
solution Eq. 7 is valid for the whole range of 
interest from if=l to -oo. 

2.3.2 Asymptotic expansion solution of the forced 
response . The asymptotic expansion solution obtain- 
ed when k = 0 was in a simple form of a divergent 
series in inverse powers of T . Because of its 
simple form it was easy to obtain the solution to 
any prescribed power, and to obtain error bounds as 
a function of the number of terms chosen. It was 
found that use of only two terms in the expansion 
gave results with very little error, and the third 
term is down by a factor of 10“8 compared to the 
first. Of course, if too many terms are taken the 
error will start to grow again, since the expansion 
is divergent. 

When stiffness ¥ is present it is no longer a simple 
matter to obtain the expansion in general form for 
any order, but the same technique nevertheless 
applies to develop the expansion term by term, and 
to bound the error at each stage. 

The expansion is obtained by first assuming that 
the time varying coefficient is sufficiently slowly 
varying that it can be considered constant for the 
purpose of generating a particular solution. This 
solution is plugged into the left hand side of Eq. 

6. The difference between the result and the right 
hand side is the forcing function for a differential 
equation whose solution is the error in the previous 
solution, and whose differential operator is the 
same as in Eq. 6. This equation is solved by the 
same method, and the procedure is repeated indef- 
initely. When three terms are developed by this 
method one obtains 


- 
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/ax \T z zot? z _ 


Hr* I 


3io tr 4 
Ck ft*) 


7 ] (">) 


A bound on the error, obtained starting from the 
variation of parameters solution, and written in 
terms o'f the -forcing function for the next error 
equation, F(f) , is 

lA< * 1 5 

The general solution of Eq. 6 obtained as the sum 
of Eqs. 7 and 10, applied to the case k=0, and using 
only two terms in the expansion, predicted well: 
the initial amplitude of the oscillations resulting 
from applying the step change in torque to the hub, 
the decay in amplitude and the increase in freq- 
uency as the spacecraft spins-up. The nonoscil- 
latory part of the solution coming from the particu- 
larsolution predicted the initially large angles 
which converge to zero as "f gets large due to 
growth of centrifugal force on the booms. These 
predictions are simple to obtain from the equations 
here, but the multibody dynamics program simulation 
of 11 minutes of spin-up, required 30 minutes of 
computatuion time on an IBM 360-91. 

3. AXIAL THRUSTING MANEUVERS 


During station keeping maneuvers it would be natural 
to require velocity and position corrections both 
in the plane of the spin and along the spin axis. 

This section is devoted to studying the effect of 
axial thrusting on the boom dynamics, both for the 
radial booms and for the axial booms. Such thrusting 
can be destabilizing to the spin axis booms, and the 
interplay and trade-offs are established here between 
such factors as: axial boom stiffness, axial boom 
length, axial acceleration magnitude, spacecraft 
spin rate, displacement of the boom axis from the 
spin axis, and misalignment of the boom direction 
from the spin axis direction. The spacecraft is 
assumed to be in a steady spin, and in the case of 
spin axis boom dynamics, the spacecraft is assumed 
to undergo a prescribed acceleration along the spin 
axis, with the boom motion having no influence on the 
hub motion. 


3.1 R esponse of radial booms to axial thrusting 

Axial thrusting directly excites the /3 spacecraft 
mode in Fig. 1. Use of Lagranges equations gives 
the boom deflection /3 in terms of arbitrary const- 
ants Ci and C 2 as 

f3 - /3 f + C, + C z 

fir = " C^lCO/fJt + ^ca*) ( (12) 

-rip = { l> + *«»«£]/ 1 (Ms - KM') ^ 


3.2 Axial boom dynamics with boom misalignment and 
offset from spin axis 


Figure 3 indicates nomenclature needed for the anal- 
ysis of the axial boom response to axial thrust. The 
coordinates T, j.'fc are inertial. The position of 
the center of the top surface of the satellite hub 
is given as ict? 1< so that the axial acceleration 
of the hub is prescribed as c. The base of the 
axial boom is displaced a distance r 0 from the spin 
axis, and hub fixed coordinates are centered at the 
base of the boom with ?' along the direction of 
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the displacement r Q . The angular orientation of_ / 
the boom is prescribed by angles 9, in the I'', K 
plane and angle 9 2 . If the boom unstressed dir- 
ection is misaligned and not parallel to the spin 
axis, the misalignment is assumed to be small 
enough for linearization to apply, and is denoted 
b y 9,0 ■ ©«• 


As before, due to space limitations, derivations 
must be deleted. Developing Lagranges equations 
for &©,, & Si and linearizing in these quantities 
as well as in 0, o , 9 m produces 

&e + G&e + k & e = [i 0 ] T + u»Ve 0 
*>© = Ue, ^®i] T 5 e c = le JO e 20 f (.3) 


©i = &i 0 + % 6 C 1=1,2. 



K . = (u, 1 -0* +■ «; e.o) 

K z = 

U>*= | + ^c./ 

w.'=(SA)Cr,/t) 

-a 1 = 


Here JTl is the radian natural frequency of the first 
mode of the boom vibration as a cantilever beam, 
and m.g is the boom mass. Consider one of these 
equations explicitly / l(1 v 

to el (»> 00 (S') VW 

£e, - 2u 0 Se 2 + [-n. l -v 0 * - + ~ ( J) e io]^6, 


+ 5“« a (r) e « %6 ^~ 

(O 


?<(£)+ <&+%)*,» 

<v) (8) 


[*»(« + ‘J?+u*-e w ) f ] 1/2 } 

Cj = tcj,, c is ] T j i = 

C tl = (2q,w,c„+ K 2 c 12 )/ CIS,- w, 1 ) 

Cj,= C V < 2 c 12 - 2 £J „w,C, 2 )/Ck,-w 1 ' 1 ) 
c 3l — C 5 w e w 2 Cin a- K 2 C31),/ (K,— ) 

0)1 s CtSjCqj— Wj C3I^/(\S> — w 2 ^ 


C, + Wj C 3 = O 

c a + Cq = -se T + se^ 


c=0 


The last two equations evaluate the arbitrary con- 
stants Cjj, C22, C32, C42 for the case of zero 
initial conditions. S© p represents the linearized 
expression for the altered equilibrium position of 
the boom due to spin and acceleration. Certain 
special cases are of interest: 

w,, 2 = \co.±^XL 1 - |f \ J 6 0 =0 or r 6 = 0 

Se p = > of •, e„= o 

H = («*•* S = ° 


Ignoring stability due to gyroscopic coupling, 
which is of no practical significance since it dis- 
appears in the presence of energy dissipation in the 
beam, stability of the differential equation occurs 
if and only if the stiffness matrix K is positive 
definite. This requires KpO and hence the square 
bracket in Eq. 14 must be positive. For any given 
boom stiffness k and boom lenght^ , there exists 
an upper limit on the allowable spacecraft spin 
rate U a combined with the allowable axial acceler- 
ation, expressed by this condition. The same effect 
is shown by the particular solution Se p . 

4. PULSED RADIAL THRUSTING 


Term (1) represents gyroscopic coupling. The brack- 
et multiplying 60| would represent the square of 
the. vibration frequency of the boom, if there were 
no coupling with the inequation, but is still 
related to the natural frequency here. This brack- 
et contains the square of the boom's natural freq- 
uency, term (2), which is decreased by the satellite 
spin rate, term (3), and by the axial acceleration 
of the spacecraft, term (4). Term (5) can either 
decrease or increase the value in the bracket de- 
pending on the sign of the boom misalign- 
ment 6,o . Term (6) represents coupling in the 
stiffness matrix due to misalignment. The offset 
of the base of the boom from the spin axis produces 
a forcing function to the differential equation, 
term (7), as does the boom misalignment © 10 , term 
( 8 ). 


The general solution for 6,, © 4 is 

©Cfc) — O 0 + £©j» + C,AAV-W,-fc 4 C 2 Ctro. "L 

+ vl 2 t + Cg cw 


$e Pi = W-ul i <£ + £ e l0 ) - e* 

W-<-£+^e, 0 f- 

se p = [ 1 6 

l(o*-u?4coj-e 10 )*- 20 


Os) 


W l,2 0),’'- <^2 + CJj Q, 0 ) + 


Because the spacecraft is spin stabilized, station- 
keeping and trajectory modifications requiring cor- 
rections in the plane of the spin are accomplished 
by turning on a radial thruster for a fraction of 
each revolution of the hub. Here we study the 
dynamic response of the spin axis booms to this 
repetitive transverse excitation at their bases. 

As in the previous section we assume that the hub 
motion is unaffected by the boom vibrations. An 
offset r 0 of the boom from the spin axis is included 
again, but misalignment is ignored, i. e. 8 o =0 . 

Figure^ 4 describes the nomenclature for the problem. 
Axes i , , , k, are inertial coordinates, axes T 2 , 

j 2 , k 2 are centered on the top of the spacecraft 
with along the spjn axis^nd along the thrust 
direction, and axes 13, $ , Tj are rotated with 
respect to the hub coordinates by angle If about 
the spin axis and translated along the resulting T 
direction by r 0 so that the origin of the coordinates 
is at the base of the boom. The path of the origin 
(>2 of the hub coordinates executes a prescribed 
zig-zag trajectory in the i, , T| plane as a result 
of the pulsed radial thrusting. 

Developing Lagranges equations for ©, and ©z> 
linearizing, and def i m ng 6 = 0, + 1 B 2 , t=^T , 
results in the differential equation 

9 + ac%ie + (.s?-oi)@ (it) 

= -liCA.+ iAl) + I Sl 5 
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A, = (F/m) co-a !f 

A a - (F/m) if 

(O 4or ■Let*'T J MT+t] , 

F= te UTvr^WOT] 

Here A^, Ag are the components of the inertial 
acceleration of the origin 03 in 3 axes, F is the 
thrust magnitude, m is the total spacecraft mass, 
and T is the period of revolution of the hub. F 
is pulsed as indicated. 

The solution of the homogeneous equation is 

0„ = C H (in.t'l + tViexp(in a tJ (n) 

XL, = XL-<*>o 
XL a ~ ~S1~ ^0 

The process of finding a particular solution is 
somewhat complicated. The forcing function is 
piecewise constant, and within each thrust-on or 
thrust-off interval a constant particular solution 
would apply. But using such a combined solution 
is simplistic in the sense that it does not cor- 
respond to a solution of the differential equation 
for any initial conditions. The particular solut- 
ion for any chosen initial conditions would norm- 
ally require solving the differential equation on 
the first interval from 0 to t , evaluating the 
arbitrary constants, and then finding ectl, elx). 
These would be used as initial conditions for the 
arbitrary constants in the solution from t to T 
using the constant particular solution. The 
process would have to be repeated indefinitely to 
obtain the solution for all t. 


However, the forcing function is a periodic funct- 
ion, which suggests that there should be a periodic 
particular solution for properly chosen initial 
conditions. By finding these initial conditions 
one only needs to find the particular solution 
over one period, and then it is known for all t. 
Hence, we look for a solution in the following form 
where the arbitrary constants are determined by the 
given boundary conditions: 


6M> = c*exp (.ifi,*) + T\exp tt F K 3 o 4 t < X 

9(fcl - C Q exp(iiL,t) +■ (i.n.j-0 “ 0<i + 1 '0 

K j ■= (Va) Aj / [X cn.*-otf )} i ='j 2 - 
Kj =(3/0 (aJ-r./uc^-Aio)] 


0 C?") = ©Ce + ) 

©ft-) = Gd + ) 

&(0) = ©CT) 

ecol = 0CT) ('8) 


The general solution for the boom response to axial 
thrusting, after employing the particular solution 
resulting from the above and after evaluating the 
arbitrary constants in terms of ©to) , ©to), is 

ett) = e P M + C M e ln ' t +D H e i - A * t tn) 




-f.r- nT^-t-xnT+t 


tj ij... 


6 pft) : 




-n-iJ 

-tor WT+'t O-'+OT 


© A Ct) = C fl e £iX ' Ct '" T + T5 ft e initt *" T + V <3 
Cfl= _ K - +LK * { e - ££} 1 

B XL, CK,+UO -oil* CTrt )/^(x^) 1 


' H ” K XL, -XW 

. _ n + i e(o)+Ji l 6(o)-n ,\<3 


If pulsed radial thrusting is maintained long enough 
that the steady state oscillations are reached, then 
the periodic particular solution represents the boom 
response. The complete solution is needed to eval- 
uate the initial transient response, and to deter- 
mine the maximum deflections encountered in a man- 
euver. 

Figure 5 presents numerical results using these 
equations. The thrust level produced 0.4 ft/sec 
acceleration (corresponding to 10 lbs thrust), and 
thrust was kept on for 90° of rotation. The boom 
lenght was 1 - 23 ft with a natural frequency of 
n = 2.756 rad/sec (corresponding to a 1-1/8 inch 
STEM boom). The boom offset r 0 was chosen as 0.1 
ft, which was considered large but possible, 
spacecraft spin was &J„ = 20 rpm = 2.09 rad/sec, 
and the thrust angle IT = 0) The top two graphs in 
the figure give the response of 0,1-0 and Sitt) 
for zero initial conditions. The second two graphs 
give the periodic particular solution for 3, and 
6s, while the final graph gives the associated 
thrust-on intervals. Although the vertical axes 
are labeled as angles, the units on the axes are 
converted to tip deflection in feet. These results 
agree well with an independent analysis using a non- 
linear multibody dynamics program. 
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